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ABSTRACT 


This  report  shows  that  the  mean  wind  through  a  given  altitude 
layer  can  be  determined  from  balloon  tracking  data  without  computing 
wind  values  for  individual  points  within  the  altitude  layer.  The 
assumptions  involved  in  the  procedure  are  similar  to  the  assumption 
made  in  numerical  differentiating  techniques,  i.e.j  the  position  can 
be  approximated  by  a  polynomial. 

Two  such  techniques  are  investigated.  One  uses  a  linear  squstio 
to  approximate  the  position;  the  other  uses  a  cubic  equation.  These 
techniques  were  applied  to  balloon  data  from  2?  balloon  tracks  and 
compared  with  results  obtained  from  numerical  differentiation.  It 
is  shown  that  there  is  very  little  difference  aiaong  the  results  of 
the  various  calculation  procedures. 
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INTRODUCTION 


Standard  techniques  for  wind  measurements  abo^e  the  level  at 
which  a  fixed  instrument  (anemometer)  is  feasible  employ  the  use  of 
balloons.  As  the  ba’ioon  ascends,  its  position  is  tracked  by  theodo¬ 
lites,  radar,  or  CM).  The  tracking  da*”a  are  processed  to  obtain  the 
positic..  of  the  balloon  at  discrete  time  periods. 

To  determine  the  wind  velocity,  one  usually  numerically  differ¬ 
entiates  the  position  data  and  assumes  that  the  horizontal  component 
of  the  wind  is  equal  (at  least  in  magnitude)  to  the  horizontal 
component  of  the  baiioon's  velocity.  The  degree  of  sophistication 
involved  in  the  numerical  differentiation  technique  depends  upon  the 
frequency,  amount  of  data,  and  the  accuracy  desired.  Once  one  has 
obtained  wind  velocities  at  ''aricus  altitudes,  the  mean  wind  through 
a  given  altitude  layer  can  be  computed  as  a  simple  average.  This 
paper  presents  and  evaluates  a  simple  technique  for  calculating  thp 
mean  wind  through  an  altitude  layer,  which  does  not  require  the  com¬ 
putation  of  the  individual  wind  velocities.  This  procedure  is  well 
suited  for  real-time  computer  applications  wherein  the  available  com¬ 
putation  time  is  limited, 

A  real-time  meteorological  svs'em  has  been  developed  by  the 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range  ;WSMR } , 

New  Mexico,  for  support  of  unguided  rocket  firings  (Duncan  and 
Rachele,  19fc7);  this  system  uses  pilot  balloons  for  wind  measurements 
from  150  to  3000  meters.  The  balloons  are  tracked  by  either  radar 
or  a  semiautomatic  triple  theodolite  system.  The  system  is  current lv 
used  for  support  of  the  Athena  rocket  and  will  also  be  used  for  support 
of  tue  Aerobee  350  rocke  .rings. 

The  techniques  presented  here  are  used  for  real-time  system  wind 
calculations.  Therefore,  it  is  of  interest  to  know  not  only  estimates 
of  the  errors  in  the  computed  wind*  but  also  the  contribution  of 
these  errors  to  the  impact  dispersion  of  the  rocket.  Both  problems 
will  be  analyzed.  It.  must  be  emphasized  that  this  paper  considers 
only  those  errors  introduced  by  the  calculation  procedures  and  does 
not  attempt,  to  consider  errors  in  the  tracking  system. 


DISCUSSION 


Let  (X,  Y,  Z)  be  a  right-handed  orthogonal  topocentnc  system 
with  the  positive  2  axis  vertical  (i.e.,  along  the  outward  normal  to 


the  earth).  The  motion  of  the  balloon  can  le  expressed  by  the 
functional  relationship  f(t).  The  component  motions  are  expressed 
functionally  by 

x  «  X(t) 
y  -  Y(t) 
z  *  Z(t) 

Since  f(t)  represents  the  motion  of  a  physical  oodv,  it  is  continuous 
and  differentiable  almost  everywhere.  It  follows  that  X(t),  Y(tj  and 
Z(t)  are  continuous  and  differentiable  almost  everywhere.  It  will 
be  assumed  for  convenience  that  2(t)  is  a  monotone  nondecreasing 
function  of  t. 


The  computation  the  X  component  of  the  mean  wind  thr.ugh  a 
given  altitude  layer  [Z(t.),  Z(t.)]  is  equivalent  to  the  computation 

1  i. 

of  X(t)  for  t.  <  t  <  t . .  This  mean,  which  will  be  denoted  bv  W  ,  is 
l  —  —  4  x 


W 

x 


t2'tl 


t2 

/  X(t)dt 


(Xtt,) 


(1) 


Thus,  if  one  can  determine  the  functions  X(t),  Y ( t ) ,  and  Z(t),  it  is 
a  simple  matter  to  compute  the  mean  wind  speed  for  a  given  altitude 
layer.  Unfortunately  it  is  extremely  difficult,  if  not  impossible, 
to  determine  these  functions.  This  suggests  the  determination  of  an 
approximating  function.  The  standard  technique  is  to  choose  a 
polynomial  of  degree  n  as  the  approx imat i ug  function  and  determine 
the  coefficients  of  the  polynomial  uy  a  curve-fitting  technique. 

When  the  computations  are  to  ne  used  in  real-time  applications, 
it  is  desirable  to  choose  n  as  small  as  accuracy  requirements  will 
polynomials  of  degree  less  than  or  equal  to  4  will  be 
considered  herein. 


In  the  remainder  of  the  discussion  u  will  be  assumed  that  the 
observations  are  equally  spaced  timewise.  It  will  also  be  assumed 
that  :,r  -  1  and  t  j  *  -t  j.  This  latter  assumption  is  nonrestrict  ive 
since  it  results  from  a  simple  change  of  variable.  Consider  the 
polynomial  of  degree  n 

no  •”  •  Vi  - *  V  *  V 

Equation  (1)  yields 

w  «  a  for  n  *  1,  2 

X  i 


x 


tor  n  « 


3  or  4 . 


The  coefficients  in  equations  (3)  can  he  determined  by  least- 
squares  techniques.  Substitution  cf  the  values  so  obtained  into  (3), 
after  considerable  simplification,  yields 

W  -  12  £  t.x/(N(K2  -  1))  for  n  *  1,  2  (•') 

10(3N3-1SN2^N»S1)  Itx  *801 t^x 

W  a  - 1— - - - - - ♦  - X - - - 

X  N(N  -1)(N*2)(9-N4)  N(N‘-1)(N*2)(N*3)  (5) 

for  n  *  3,  4 

where  N  is  the  number  of  sample  points  in  the  layer,  and  the  summation 
tend  over  the  interval. 

The  calculations  required  for  n  =  1  ~r  2  are  quite  simple, 
especially  when  compared  to  that  required  for  the  n  -  3  or  4  case. 

EMPIRICAL  EVALUATION  O'-  THE  MEAN  WINDS 

The  errors  associated  with  equations  (4  and  ( S )  can  he  evaluaTeJ 
theoretical!)  by  assuming  various  functional  forms  for  the  bai leer 
path.  However,  the  magnitude  cf  the  number  of  conceivable  paths 
suggests  that  empirical  evaluations  may  be  more  feasible,  1'oward 
this  end,  2"  100-gram  pilot  balloon  track*  were  evaluated.  The  >•  olloo 
were  tracked  by  a  Contraves  c ’ netheodol i te  system.  The  reduced  data 
consisted  of  position  and  velocity  d-’ta  computed  at  the  rate  of  one 
sample  point  per  second.  The  velocity  component  s  were  c or pu t cd  b\ 
a  numerical  differentiation  technique  described  bv  Comstock  et  a  1 . 
(1964) ,  The  average  velocities  .denoted  'Jv}  were  used  as  a 
basis  for  the  comparisons. 

The  27  balloon  soundings  were  separated  into  two  groups .  In  the 

first  group,  consisting  of  13  soundings,  equations  ;  4 :  «,nd  •  ■  and 

their  analogs  for  the  y  component  were  evaluated  tor  adjacent  height 

layers  of  150  meters.  Adjacent  height  layers  of  60  meters  were  used 

for  the  data  in  the  second  group.  The  results  for  the  x  •o-.ponent 

are  denoted  K  and  H  respectively.  The  wind  profiles  resulting 

x .  x , 

i 

from  the  three  calculation  procedures,  for  a  typical  balloon  trac>, 
are  shown  in  Figure  1.  It  is  easv  to  see  that  the  differences 
between  any  two  of  these  wind  profiles  arc  quite  xr.a 1 1 . 


APPLICATION  TO  ROCKET  IMPACT  PREDICTION 


sti  •’ 


$  t 

I  i 


The  differences  AXj  -  -  Ux,  Ax2  -  -  Ux»  and  Axs  -  Wx  - 

Wx2  an<*  *^eir  an»logs  for  the  y-component  were  computed  for  each 
point.  These  data  were  used  to  compute  RMS  values  o.  ,  a,  .a  . 

AXj  AX2  AXj 

°Ay  *  °Ay  *  an<*  °Ay  *0T  eac^  h*ll°0n  track;  ensemble  values  were 
1  ^  3 

computed  for  the  group  I  data  and  the  group  II  data.  Those  results  are 
listed  in  Table  I. 


Since  the  winds  computed  by  equations  (4)  or  (5)  are  used  for ~ 
rocket  impact  prediction,  it  is  natural  to  investigate  the  contribution 
of  these  errors  to  the  error  in  the  predicted  rocket  impact.  Although 
this  question  cannot  be  answered  in  general,  it  can  be  investigated 
for  specific  rocket  configurations.  Such  investigations  should  provide 
insight  toward  the  answer  for  other  rockets  of  similar  ballistic 
characteristics. 

The  rocket  impact  dispersion  (change  in  impact  from  nominal), 
for  a  given  wind  profile,  can  be  computed  by  rocket  trajectory  simula¬ 
tion  [Cochran  et  al.,  1966],  The  rocket  impact  dispersion  for  each 
of  the  wind  profiles  discussed  in  the  preceding  section  was  computed. 
Two  different  rockets  were  considered:  the  Aerobee  150,  a  high- 
altitude  research  rocket;  and  the  Athena,  a  reentry  research  rocket. 
These  results  are  listed  in  Tables  II  and  III. 


CONCLUSIONS 


Techniques  for  computing  wind  velocities  from  pilot  b&xloon  data 
'  have  been  presented  and  compared.  Three  methods  were  discussed:  (1) 

a  least-squares  linear  fit  to  the  position  data;  (2)  a  least-squares 
j  cubic  fit;  and  (3)  averaging  of  point  by  point  results  of  numerical 

differentiation. 

^he  three  techniques  were  compared  by  applying  each  computation 
to  27  different  balloon  tracks.  The  results  for  one  such  track,  which 

t  , 


5 


is  considered  typical,  are  shown  in  Figure  1.  It  is  easy  to  see  that 
there  is  very  little  difference  among  the  three  wind  profiles.  The 
RMS  values  presented  in  Table  I  also  indicate  at  the  three  compu¬ 
tational  proce  res  yield  small  differences  among  the  resulting  wind 
profiles. 

The  various  wind  profiles  were  used  to  compute  the  impact  disper¬ 
sion  of  the  Aerobee  and  Athena  rockets.  The  results  are  .hewn  in 
Tables!!  and  III.  The  differences  among  the  three  wind  profiles 
are  considered  to  be  quite  small,  especially  when  compared  with  such 
contributors  to  dispersion  as  thrust  misalignment,  wind  variability, 
etc . 

The  authors  contend  that  the  differences  among  the  tore-  compu¬ 
tational  procedures  are  insignificant.  This  is  especially  true  if 
the  computed  wind  data  are  to  be  used  for  rocket  impact  prediction. 

Any  of  the  three  technique,.  be  applied;  i t owever” sfTrprfl rtTV~17?~The 
calculations  suggests  the  line*,  fit,  and  this  method  is  recommended 
for  use  in  real-time  rocket  impact  prediction. 


TABU'.  I: 


RMS  Differences, Meters  Per  Sec, 


i<un  No. 

*  V 
‘"*'1 

~  ix., 

5» 

'  7,x  , 

j. 

“  7-v . 

-  i 

"  6y, 

y  7, 

I 

1 

* 

0.18 

n  i 

\  19 

0.16 

0,24 

0,17 

I 

\ 

<V 

0.11 

0.10 

0.1? 

0.11 

0.07 

0,13 

I 

•7 

o 

0. 14 

0.1S 

0.19 

0.23 

0.10 

0.21 

* 

« 

<4 

0,1? 

0.11 

0.18 

c.r 

0  i 

0.19 

T 

5 

0.1S 

0.08 

0.19 

i  -• 

V  »  i  C. 

0.10 

0.14 

I 

6 

0.14 

0.10 

0.19 

0,20 

0.16 

0,20 

T 

i 

7 

0.20 

0.13 

0.24 

0.19 

g  l  7 

0.16 

I 

s 

0,25 

0.13 

0.25 

0.33 

0.17 

0.42 

r 

i 

9 

0.09 

0,08- 

0.12 

0,14 

0.07 

0.16 

~r~ 

TU 

..  jj-jg- 

0,23 

o.ir 

0.25 

0,18 

0.27 

i 

11 

0.16 

0.13 

0.22 

0,11 

0.09 

0.14 

r 

i 

12 

0.10 

0.05 

0.14 

0.11 

0.07 

0.09 

I 

13 

0.11 

0.07 

0.1? 

0.1? 

0.10 

0.14 

I 

14 

G.I9 

0.18 

0.14 

0.07 

0.09 

0.10 

I 

IS 

0.1? 

J.  10 

0,20 

0,14 

0.07 

0.1S 

I* 

o,*s 

0.11 

0.18 

0.16 

0.10 

0.17 

II 

1 

0,14 

o.os 

0.16 

0.10 

0.06 

0.12 

II 

7 

0.28 

0.24 

C.31 

0,36 

0.31 

0.30 

II 

3 

0.S2 

0.39 

0.51 

0.40 

0.3S 

0.24 

II 

4 

0.33 

0.20 

0.32 

0.42 

0.33 

0.30 

II 

S 

0.23 

0. 23 

0.19 

0.20 

0,22 

0,15 

II 

6 

0.21 

0.  is 

0.16 

0,16 

0.13 

0.13 

II 

7 

0 . 35 

0.25 

0,2c 

0,24 

0.14 

0.21 

II 

8 

0,16 

0.13 

0,15 

0,13 

0.14 

0.14 

II 

9 

0.2? 

0,17 

0.22 

0,26 

0.22 

0.25 

II 

10 

0.13 

0.06 

0.18 

0.09 

0.06 

0.13 

II 

II 

0.23 

0.18 

0.29 

0.26 

0.29 

0.26 

11 

12 

0.2? 

0,21 

0.21 

0.25 

0.17 

0.20 

II 

* 

0.28 

0.22 

0.23 

0.26 

0.23 

0.22 

*  Ensemble  Estimates, 
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I:  AFROBEE  IMPACT  DI5"  RSION  FOR 


X-Component 

Meters 


Profile 

Prof i le 

Prof  i 

W 

K 

If 

*1 

X2 

X 

14064, 

14236. 

14593 

35761  . 

36486. 

35282 

31  i  1 6. 

31431. 

31350 

33439. 

33509. 

33373 

32618. 

32S82 , 

32660 

-8136. 

-7601 . 

-7766 

-3983. 

-3501. 

-3726 

-5282. 

-4715. 

-4493 

7081. 

7060. 

6962 

-6012. 

-5895. 

-5887 

21525, 

2160S . 

?  2320 

27367. 

27119, 

27208 

24SSS. 

2561 9, 

z5366 

18056. 

18370. 

17933 

24870. 

24094. 

24326 

5783. 

4921 . 

5482 

••.'•240. 

-2933. 

-3125 

-8C73. 

-7457. 

-7973 

-16415. 

15662 . 

-15304 

12569. 

14533. 

14873 

-1911. 

-2004. 

-2009 

16285. 

19460. 

21791 

29869. 

30494. 

30871 

224S1 . 

23014. 

24752 

16810. 

17568. 

16908 

3  2°%. 

30S68 . 

31087 

56423. 

57807. 

58645 

DIFFERENT  WIND  PROFITS 


Y-Comportent 

Meters 


Prof i 3  e 

Profile 

Prof  x 

W 

N 

U 

yl 

4. 

y 

4851. 

5145. 

5275. 

33174. 

33831. 

33007, 

27778. 

28098. 

28122. 

33210. 

32624. 

32709. 

31592. 

31425. 

31259. 

51649. 

..1847, 

52465. 

50845. 

S0958 . 

50097. 

50586. 

52291. 

51424, 

46872. 

47202. 

47072. 

47710. 

47683. 

47072. 

-2742, 

-2235. 

-2436. 

-1675. 

-1340. 

-1283, 

-3052. 

-2769. 

-2265. 

-7516. 

-8036. 

-7549. 

-5772. 

-  5 / 53 . 

-5759. 

1241, 

930. 

1350. 

80445. 

80877. 

81326. 

74906, 

75323. 

74445. 

S237S „ 

81500. 

825S6. 

68361. 

68203. 

6777'. 

72171. 

71559. 

7134b. 

-7429. 

-5936. 

-<4310. 

875. 

854. 

1059. 

-663. 

-3526. 

-3007. 

-7988. 

-6136. 

-5210. 

-3025. 

-3213. 

-6156. 

35628. 

35165. 

34270. 

TABLE  III:  ATHENA  IMPACT  DISPERSION  FOR  DIFFERENT  WIND  PROFILES 
X -COMPONENT  Y- COMPONENT 


Meters  Meters 


Run  No, 

Profile 

Pro4'4  le 

Profile 

Profile 

Profile 

Profile 

W 

W 

U 

W 

W 

U 

X1 

%2 

X 

*1 

^2 

y 

I  I 

2S391 « 

28973. 

3005S. 

-540'* , 

-5213. 

-4913. 

I  2 

126142. 

127792. 

126304. 

14580. 

14509. 

14305. 

I  3 

121369. 

122123, 

121209. 

13358, 

13449. 

13434. 

I  4 

121030. 

120757. 

120077. 

1SS89, 

ioOS9. 

15327. 

I  S 

126823, 

126182. 

126308, 

15541. 

15566. 

15503. 

I  6 

-23840. 

-22386. 

-22769. 

33047. 

33080, 

32914. 

I  7 

-14813. 

«]  ..819, 

-13376. 

33780. 

33864. 

33987. 

l  8 

-13842. 

-12849, 

-32849, 

34850. 

34787. 

34787. 

I  9 

28215. 

28501. 

2*T772„ 

33496. 

33489. 

33501 . 

I  10 

-18920. 

-18251 . 

-18229, 

33659. 

33117. 

33658. 

i  i: 

64929. 

65475. 

66899, 

-10597, 

-9483. 

-9920. 

I  12 

79920. 

79194. 

79645. 

-12164. 

-11418. 

-1  lost) , 

I  13 

66594. 

71107. 

70173, 

-12800. 

-12589. 

-1,1922. 

I  14 

5133^ 

51019. 

50027, 

-  -14836. 

-14842. 

-14280. 

I  15 

67613, 

6S246. 

6S985 . 

-11759, 

-11 654. 

-13676. 

II  1 

28104. 

26432. 

2^518, 

-5772. 

-5901 , 

-5504, 

II  2 

-23608. 

-22233. 

-’1981  . 

31  ''58. 

31828. 

31797. 

II  3 

-8901. 

-9312. 

-11172, 

33680. 

33677. 

33656. 

II  4 

-16417. 

-15764. 

-16150. 

33ul9, 

33649 

33540. 

II  5 

25353. 

28095. 

29691 . 

33600. 

33474, 

33488. 

11  6 

-16931. 

-17374. 

-17981  . 

33603. 

33588. 

33623 . 

II  7 

60775. 

65983. 

6’1  69. 

-9142. 

-9647. 

-8979. 

II  S 

80881. 

81534 , 

81828, 

-10575. 

-10893. 

-10790. 

II  9 

63144, 

65 97 7. 

68L.S. 

-8959. 

-11350. 

-11522. 

II  10 

45815. 

46180. 

45013. 

-15129. 

-  14535. 

-13984. 

II  11 

71199. 

6 7 314. 

67214. 

-11277. 

-1225- . 

-12665. 

II  12 

156728. 

156335. 

156934. 

9673. 

96  1 . 

9809. 
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This  report  shows  that  the  mean  wind  through  a  given  altitude  layer 
can  be  determined  from  balloon  tracking  data  without  computing  wind  values 
for  individual  points  within  the  altitude  layer.  The  assumptions  involved 
in  the  procedure  ar  similar  to  the  assumption  made  in  maaerical  differ¬ 
entiating  techniques,  i,e.,  the  position  can  be  approximated  by  a  polynomial. 
Two  such  techniques  are  investigated.  One  uses  a  linear  equation  to 
approximate  the  position;  the  other  uses  a  cubic  equation.  These  techniques 
were  applied  to  balloon  data  from  27  balloon  tracks  end  compared  with 
results  obtained  from  numerical  differentiation.  It  is  shown  that  there 
is  very  little  difference  among  the  results  of  the  various  calculation 
procedures. 
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